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ABSTRACT 



An overview and investigation of the more popular digital filter design tech- 
niques are presented, with the intent of providing the filter design engineer a com- 
plete and concise source of information. Advantages and disadvantages of the 
various techniques are discussed, and extensive design examples used to illustrate 
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I. INTRODUCTION 



Digital filter design today is a rapidly growing field with myriad applications 
in the areas of signal processing, image processing, filtering, prediction, and estima- 
tion. Extensive research has resulted in a plethora of available information, which 
can be overwhelming to the engineer presented with a specific design problem. As 
the title suggests, the purpose of this thesis is to present an overview of the more 
popular design techniques in an attempt to consolidate the information available in 
the literature, and provide an easily readable reference manual. Detailed examples 
are given to illustrate the application of selected techniques to specific filter design 
problems. 

Due to the wealth of information available, a summary and investigation of 
every filter design technique is not possible. As stated the more popular meth- 
ods are outlined in detail, however, references are included as sources of further 
information. 

The survey is threefold, consisting first of recursive HR design methods, fol- 
lowed by nonrecursive FIR design, and finally computer-aided design (CAD). A 
summary of the chapter contents follows to enable the reader to immediately lo- 
cate the section applicable to his/her particular design problem. 

Chapter II, Recursive Filter Design, presents Butterworth, Chebyshev and 
elliptic filter design from both a traditional approach, wherein analog prototype 
filters based on design specifications are converted to digital versions using the 
bilinear transformation, and a direct design approach that eliminates the need to 
determine an analog prototype. 
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Chapter III, Nonrecursive Filter Design, includes an overview of Fourier co- 
efficient design, windowing, and a derivation and illustration of the method of 
frequency sampling. Among the design examples presented is the design of a band- 
pass differentiator and a bandpass integrator. 

In Chapter IV, Computer-Aided Design, the Optimum FIR Filter Design 
Method of Parks and McClellan (that employs the Remez Exchange Algorithm) 
is investigated and applied to a bandpass filter design problem. Also included is 
a short discussion of a method to enhance the application of this program to the 
design of high-order filters. 

To complete the chapter, the Minimum p-Error Design Method for HR Filter 
Design, that utilizes the Fletcher- Powell algorithm, is presented. Furthermore, a 
discussion of the advantages and disadvantages of these two iterative techniques is 
included. 
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II. RECURSIVE FILTER DESIGN 



A. INTRODUCTION 

The recursive realization of digital filters is advantageous, in that the desired 
frequency response can be obtained using a lower order filter than if a nonrecursive 
realization were used (assuming linear phase is not required). This is because the 
filter frequency response is influenced by both the poles and zeros of the filter 
transfer function, whereas the frequency response of a nonrecursive realization is 
determined only by the filter’s zeros. 

Traditionally, the design of recursive digital filters involves the determination 
of an analog prototype using one of several methods: Butterworth, Chebyshev or 
elliptic, to name a few, and converting these to a digital version using impulse- 
invariant design or bilinear transformation methods, the latter being the more 
popular. 

This chapter first presents a review of the traditional design techniques us- 
ing examples involving Butterworth, Chebyshev and elliptic filters; then a Direct 
Design method is presented, whereby the requirement for determining an analog 
prototype is eliminated, thus reducing the algebraic complexity of recursive filter 
design. 

B. TRADITIONAL DESIGN TECHNIQUES 

Recursive filter design techniques, as stated in the introduction, involve deter- 
mining an analog prototype filter from given filter specifications, and then convert- 
ing the prototype to a digital version using a bilinear transformation. Since this is 
a well known design procedure [1], details of the derivation will not be considered. 
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However, as a quick review, three examples involving the design of Butterworth, 
Chebyshev and elliptic filters will be given. 

It should be noted that the tables of analog prototype transfer functions found 
in this section (Tables 2.2, 2.3 and 2.4), all have a critical frequency of unity 
( w c = 1). Thus, to obtain a filter transfer function with a different critical fre- 
quency based on these prototypes requires the use of an appropriate analog fre- 
quency transformation (Table 2.1). 

For example, suppose a lowpass analog filter with a critical frequency of w c 
is desired. The following relationship applies: 



HlpU w c) = Hip p {j\) 



( 2 . 1 ) 



where Hi Pp is the prototype filter with a critical frequency of w c = 1. 

To convert the prototype transfer function to the desired transfer function, 
^Lp(«s)? the s in the lowpass prototype filter is replaced by s/w c as: 



H LP (s) = H LPp {s) 



( 2 . 2 ) 



Table 2.1 (Reference 1) lists the appropriate frequency transformations for high- 
pass, bandpass and bandstop filters. A summary of the design procedure follows. 

1. Bilinear Transformation Design 

From given analog filter specifications, determine the approppriate s- 
domain design technique (i.e., Butterworth, Chebyshev or elliptic). 

Design the analog prototype filter in accordance with the technique se- 
lected above. This involves the translation of the filter specifications to those of 
a lowpass prototype. The lowpass prototype is then designed according to the 
translated specifications. 
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Transform the lowpass prototype to the desired lowpass, highpass, band- 
pass, or bandstop filter according to Table 2.1. 

Transform the analog filter transfer function to a digital version using the 
bilinear transformation. 



H(z) = H(s) 



(2.3) 



Example 2.1 Butterworth Filter 

Design a digital filter that is flat in the passband from 0 to the 3 dB cutoff 
frequency, f c , of 2 kHz. For frequencies greater than 4 kHz, the attenuation should 
be at least 10 dB. The sampling r*\tc is 20 kHz. 

Step 1 : A Butterworth design is called for because a flat passband is desired. 

Step 2 : Convert the critical analog design frequencies to digital. 



Sampling frequency: f s = 20 kHz ==> w s = 40 x 10 3 7 t rad /s 

Cutoff frequency: f c = 2 kHz ==> w c = 4 X 10 3 7r rad/s 

Stopband frequency: f a = 4 kHz ==> w a = 8 x 10 3 7r rad /s 

= w cT = W c /f, = ( ^ 0 X x 1 f 0 ) 3 T = 0.2* rad 

e 'a = w aT = w a /f, = ^ 0 >< x 1 ° 0 3 ,r = 0.4* rad 
Step 3 : Prewarp the analog frequencies to yield the desired digital frequencies. 

w' c = tan (0^/2) = tan(0.1;r) = 0.325 rad/s 
w' a = tan(^/2) = tan(0.27r) = 0.726 rad/s 



Note: Prewarping is done to ensure the analog filter will ultimately 

yield a digital filter with the correct critical frequencies. This is best 
visualized in the following graph (Figure 2.1) of the relationship between 
the analog filter frequencies and the desired digital frequencies, 
i.e., w' n = tan (6 n /2). 
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TABLE 2.1 



ANALOG FILTER FREQUENCY TR ANSFORMATIONS 
(after Ref. [1]) 



Filter Type 



inform s-form 

to find (to find H{s) from H lPf (s) 

replace s in prototype 
with) 
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Normalizing to the lowpass prototype, according to Table 2.1: 



w = 0.325 — ► w =1 rad 

c c p 

w' a = 0.726 — * w ap = 0.726/0.325 = 2.234 rad 




Figure 2.1. Relationship Between Analog and 
Digital Filter Frequencies 



Note: Normalizing is done to transform the prewarped critical frequencies 
of the analog filter to their relative equivalents in a lowpass prototype filter with 
a critical frequency of w Cp = 1. This is done to enable the designer to take 
advantage of previously compiled tables or frequency response curves corresponding 
to the transfer functions of these lowpass prototypes. Thus, a prototype filter can 
be selected whose frequency response characteristics have the same shape as the 
filter that is being designed. The transfer function for the selected prototype is 
then converted to a transfer function that exhibits the desired frequency response 
characteristics using an appropriate substitution (Step 5). Figure 2.2, depicts the 
process of normalizing. 
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Step 4 : 



Step 5 : 



Step 6 : 




Find the order N of the lowpass Butterworth filter prototype. The order 
N to provide a gain of A/dB at an angular frequency w a is 



N = 



log 10 (10 ( ~ MdB/10) - 1) 
21og ]0 w a 
Iog, 0 (10‘ - 1) 



; A/dB = — lOdB 



(2.4) 



21og 10 2.234 

— 1.37 (Since N has to be an integer choose, N =2) 

From Table 2.2 determine the transfer function for the lowpass filter. 



Hlp{- s) = Hlp p (s] 

1 



= H LPp (s) 



s“ 4- V2s + 1 
0.106 



(.2-5) 



s 2 + 0.460s + 0.106 

Determine the transfer function for the digital version of the filter using 
the bilinear transformation. 



Hlp(z) = Hlp(s) 

z 2 + 2 : + 1 



0.06S(r + l) 2 



( 2 . 6 ) 



14.S4z 2 - 17.0 Z + 6.14 z 2 - 1.14z + 0.413 
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TABLE 2.2 



BUTTERWORTH PROTOTYPE COEFFICIENTS 

(Table after to Ref. 4) 



N 


a i 


a 2 




a 4 


as 


a 6 


a 7 


<*8 


1 


1.0000 
















2 


1.4141 


1.0000 














3 


2.0000 


2.0000 


1.0000 












4 


2.6131 


3.4142 


2.6131 


1.0000 










5 


3.2361 


5.2361 


5.2361 


3.2361 


1.0000 








6 


3.8637 


7.4641 


9.1416 


7.4641 


3.8637 


1.0000 






7 


4.4940 


10.0978 


14.5918 


14.5918 


10.0978 


4.4940 


1.0000 




8 


5.1258 


13.1371 


21.8462 


25.6884 


21.8462 


13.1371 


5.1258 


1.0000 



LPp ^ ^ 1 -f ajs -f a2«s 2 + • . . -f a^s N 
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Step 7 : Obtain a plot of the digital filter frequency response to see if the design 

specifications have been met. 

Figure 2.3 shows that the design does indeed comply with the given filter 
specifications, in that the cutoff frequency is 0.628 = 0.27T rad and the stopband 
frequency is 1.257=0.47r rad with gains of —3 dB and —14 dB, respectively. 
Example 2.2 Chebyshev Filter 

Design a digital bandpass filter with the following specifications: 

• 1 dB ripple in the frequency band 600 to 900 Hz, 

• sampling frequency of 3 kHz, 

• maximum gain of —40 dB for 0 < / < 200 Hz. 

Step 1 : Convert the critical analog design frequencies, w,-, to the corresponding 

digital frequencies, 

Sampling frequency: f 9 = 3 kHz 

Lower ripple passband frequency: ft = 600 Hz =>» wt = 1200?r rad / s 
Upper ripple passband frequency: f u = 900 Hz => w u = 18007T rad / s 
Stopband frequency: f a = 200 Hz => w a = 4007T rad /s 

9' ( = w e T = w e / f, - ^qqq^ = 0 An = 1.26 rad 
6' u = w u T = w u /f, = = 0.6tt = 1.88 rad 

9' a = w a T = Wa/fs = = 0.133tt = 1.54 rad 

Ripple band center frequency : 0q = = 0.49?r = 0.418 rad 
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Step 2: P re warp the analog frequencies so that the desired digital frequency 



characteristics will be achieved. 

w' ( = tan(0j/2) = tan(1.26/2) = 0.729 rad /s 

w' u = tan(^/2) = tan(1.88/2) = 1.369 rad/s 



w' a = tan (0^/2) = tan(0.418/2) = 0.212 rad/s 

w' 0 = \J = \/(0.729)(1.369) = 1.00 rad /s 
Ripple band is: B' = w' u — w' ( = 1.369 — 0.729 = 0.64 rad/s 

From Table 2.1 convert the bandpass filter frequencies to the correspond- 
ing lowpass prototype filter frequencies. 



WLPp = 



W BP - W 0 2 
B'iv'bp 

w bp ~ (RQQ) 2 
(0.64)^ 



(2.7) 



BP 



w 



w' 



Step 3 : 



BP w LPp 

0.729 -1 

1.369 1 

w' a 0.212 -7.04 

w' 0 1.000 0 

Determine the order of the Chebyshev prototype filter that meets the 
design requirements. 

For a lowpass Chebyshev filter the magnitude - squared characteristic is 



given by: 



H LPp (jw) 



( 2 . 8 ) 



1 + e2 C%(w) 

where Cn(w) is an N tk - order Chebyshev polynomial and e indicates the degree 
of ripple. 

In this example the —40 dB gain translates to M = 10“ 2 for the desired 
filter frequency response in the stopband, i.e., we want < 10~ 2 or 
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|iJ(ju;)| 2 < 10 -4 . The design specification for 1 dB ripple in the ripple band 
corresponds to a value for e of 0.2589. 

Making an initial guess for the filter order of N = 2 yields: 



2 



H 2 (jw) 



1 

1 + 0.2589 (2wl -l) 2 

1 

1 +0.2589(99.12- l) 2 



w a = 7.04 

- = 4.01 x 10~ 4 



(2.9) 



which is not less than 10 4 . 

Proceeding further with N = 3 yields: 

1 






1 + 0.2589 (4iy^ — 3w a ) 
1 



1 + 0.2589(1395.65 — 21.12) 



u> a =7 .04 

= 2.044 x 10" 



( 2 . 10 ) 



which is less than 10~ 4 . Therefore, the design needs can be met using a third-order 
Chebyshev lowpass prototype filter. 

Step 4 : Obtain the transfer function of the third order lowpass prototype filter 
with 1 dB ripple. Looking at Table 2.3, yields the following prototype 
transfer function: 



Hlp p (s) 



0.491 

s 3 + 0.988s 2 4- 1.238s + 0.491 



( 2 . 11 ) 



Since this is an odd order filter, the constant 0.491 in the numerator was 
selected to make |if(j0)| = 1. Recall that for Chebyshev filter design, the constant 
K in the numerator of the lowpass prototype filter is selected on the following basis: 



for N odd; 
for N even; 



H(j 0) 

H(J 0 ) 



= 1 



[1 + e : 



211/2 



( 2 . 12 ) 
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TABLE 2.3 

CHEBYSHEV - PROTOTYPE DENOMINATOR POLYNOMIALS 
(Table due to Reference 4) 

0.5-dB Ripple Chebvyshev Filter (e = 0.3493) 



n 

1 *4-2.863 

2 ** + 1.425* + 1.516 

3 * 3 + 1.253*2 4- 1.535* + 0.716 - (* + 0.626X** 4- 0.62 6* + 1.142) 

4 j4 + l - 197,J + 1-717*2 4- 1.025* 4- 0.379 - (** + 0.351* 4- 1.064X* 2 4- 0.845* + 0.356) 

5 * J - 1.17251** -v 1.9374* 3 4- 1.3096*2 4- 0.7525* 4- 0.1789 

(* - 0.3623)[(* -r 0.1120)2 -l 1.01 16 2 ][(* -r 0.2931)2 4. 0.62522) 

6 *« - 1 . 1 592* 5 - 2.I718* 4 - 1.5898* 3 + 1.1719*2 + 0.4324* 4- 0.0948 

[(* 4- 0.0777)2 _ 1 .0085 2 ][(j - 0.2121)2 4- 0.73822][(* 4r 0.2898)2 4- 0.2702 2 ] 

7 *7 1.1512*6 4. 2. 4126* 5 1.8694** - 1.6479* 3 - 0.7556*2 4. 0.2821* + 0.0447 

(* - 0.2562)[(* - 0.0570)2 -t- 1.0064=][(* 4- 0.1597)2 4- 0.8001 =]{(* nr 0.2308)2 4- 
0.44792] 

8 *« - 1.146I* 7 -4 2.6567*6 4- 2.1492* 3 4- 2.1840** 4- 1.1486*3 4- 0.5736*2 + 

0.1525* - 0.0237 

[(* -r 0.0436)2 f 1.0050 : ][(* -r 0.1242)2 -r 0.8520 3 ][(* 4- 0.1859)2 -r 0.5693 2 ] 

[(* -r 0.2193)2 - 0.19992] 

9 *’ 4- 1.1426** 4- 2.9027* 7 4- 2.4293** 4- 2.7815*3 4- 1.61 14** 4- 0.9836* 3 4- 

0.3408*2 4- 0.0941* - 0.0112 

(* - 0.1984 )[(* -T 0.0345)2 - 1.00402][(* 4- 0.0992)2 + 0.88292][(* 4- 0.1 520) 2 4- 
0.6553 2 ][(* - 0.1864)2 ^ 0.34872] 

10 *»o - 1.1401*’ -p 3.1499** 4- 2.7097* 7 -4 3.4409** -r 2.1442* 3 4- 1.5274** 4- 

0.6270* 3 t- 0.2373*2 j- 0.0493* - 0.0059 

[(* - 0.0279)2 ~ 1 .0033 : ][(* - 0.0810)2 _ 0.905P][(j + 0.1261)2 4. 0.71832] 

[(* - 0.1589)2 _ 0.46 1 2-][(* 4- 0.1761)2 4- 0.15892] 



25 



TABLE 2.3 (Continued) 

CHEBYSHEV - PROTOTYPE DENOMINATOR POLYNOMIALS 
„ 1.0-dB Ripple Chebyshev Filter (e = 0.5089) 

1 * + 1.308 

2 *= - 0.804* - 0.637 

3 J J - 0.738* 2 - 1.022* - 0.327 = (* - 0.402N* 2 - 0.369* - 0.886) 

4 j 4 - 0.7 1 6* 3 — 1 ,256s ‘ — 0.517* - 0.206 - (*= — 0.210* - 0.928)(* = — 0.506* - 0.221) 

5 * 3 - 0.7065* 4 - 1 .4995* 3 - 0.6935* 2 - 0.4593* - 0.0817 

(s - 0.21 83)I< j - 0.0675)= - 0.9735 =][(* - 0.1766)= - 0.6016=] 

6 J 6 - 0.701 2* 3 - I.7459* 4 - 0.S670* 3 - 0.7~15*= - 0.2103* - 0.0514 

[f* - 0.0470)- - 0.98 1 7 -][( j - 0.1283)- - 0.71S7=][(* - 0.1753)= - 0.2630=] 

7 j 7 - 0.6979* 6 - 1.9935* 3 - 1.0392* 4 - 1.1 444* 3 - 0.3S25* = - 0.16661* - 0.0204 

(s - 0.1553)[(* - 0.0346)- - 0.9867=](( * - 0.0968)- - 0.7912-][(j - 0.1399)- - 
0.4391 -] 

8 *• ~ 0.696 Is 7 - 2.2423* 6 - 1.21 17* 3 - 1.5796* 4 - 0.5982* 3 - 0.3587* = - 

0.0729* - 0.0129 

[<* - 0.0265 ) z - 0.9898=] ((* - 0.0754)= - 0.8391 =][(j - 0.1129)- - 0.5607=] 

[(* - 0.1332)= - 0. 1969=) 

9 s 9 - 0.6947*' J- 2.4913* 7 - 1.3837*' - 2.0767* 3 - 0.8569* 4 - 0.6445* 3 - 

0.16644*- — 0.0544* - 0.0051 

(* - 0.1206)[(* - 0.0209)- - 0.9919=][{* - 0.0603)= - 0.8723=]I(* - 0.0924) 2 
- 0.6474-1K* - 0.1 134 >= - 0.3445=] 

10 *'° - 0.6937*’ - 2.7406*' - 1.5557 s 7 - 2.636 3*' - 1.1585* 3 - I.0389* 4 + 

0.317B* 3 - 0.1440* 2 - 0.0233* -• 0.0032 

[(* - 0.0170)2 - 0.9935=]|(* - 0.0767)2 . 0.71 13=][(j 0.0493)= - 0.8962=] 

[(* - 0.0967)2 - 0.4567 = ][(* - 0.1072) 2 - 0.1574=] 



2-dB Ripple Chebyshev Filter (e = 0.7648) 

1 *+ 1.965 

2 *2 - 1.098* - 1.103 

3 * J - 0.988*2 - 1.238* - 0.491 = (*-*- 0.494X* 2 - 0.490* - 0.994) 

4 j 4 - 0.953* 3 4- 1.454*: -r 0.743* - 0.276 = (* = - 0.279* - 0.987)(* 2 - 0.674* - 0.279) 

5 j 3 - 0.9368* 4 - 1.6888* 3 - 0.9744* 2 j- 0.5805* 0.1223 

(* - 0.2S95)[{* - 0.0895)2 -a- 0.9901 2 ][(* 0.2342) 2 - 0.61 19 2 ] 

6 j 4 -J- 0.9282* 3 -r 1.9308* 4 - 1 .202 1 * ^ - 0.9393* 2 - 0.3071* - 0.0689 

[(* - 0.0622)2 - 0.9934=][(* 0.1699)2 - 0.7272=][(* - 0.2321) 2 - 0.2662 2 ] 

* 7 - 0.9231* 6 - 2.1761 j 3 - 1.4288* 4 - 1.3575* 3 - 0.5486*= - 0.2137* - 0.0307 
<* - 0.2054 )[(* - 0.0457)2 - 0.9953=][(* -r 0.1281)= - 0.7982 2 )[(j ^ 0.1851)= - 
0.4429=) 

8 *» - 0.9198* 7 - 2.4230* 6 - I.6552* 3 - 1.8369* 4 - 0.8468* 3 - 0.4478*= + 

0.1073* - 0.0172 

K* - 0.0350)2 - 0.9965 =]((* - 0.0997)= - 0.8448 = )[(* -u 0.1492)= - 0.5644 2 ] 

((* - 0.1759)= - 0.1982=] 

9 * 9 - 0.9175*' - 2.6709* 7 - 1.8813** - 2.3781* 3 - 1.2016* 4 - 0.7863* 3 j- 

0.2442*2 - 0.0706* - 0.0077 

(* - 0. 1593)[(* - 0.0277)2 - 0.9972 2 )[(j - 0 0797)= -u 0.8769=)I(* - 0.1221)= - 
0.6509=]f(j _ 0.1497)= - 0.3463=) 

10 *io _ 0.9159*’ - 2.9195*' - 2.1079* 7 - 2.9815* 4 - 1.6830* 3 - 1.2445* 4 - 
0.4554* 2 - 0.1825* 2 - 0.0345* - 0.0043 

K* - 0.0224)2 - 0.9978 2 ][{* * 0.I013) 2 - 0.7143=)((* - 0.0651)= - 0.9001=) 

[(* - 0.1277)2 - 0.4586=)[(* - 0.1415)2 - 0.1580=] 
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